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Introduction 

The  main  achievements  of  the  activity  carried  out  under  this  grant  have  been 
an  understanding  of  the  mechanism  of  rain  noise  and  of  the  role  of  bubble 
clouds  in  ambient  noise  generation  at  low  frequencies. 

The  first  problem  has  been  investigated  theoretically  in  conjunction  with 
experiments  carried  out  at  the  University  of  Mississippi  by  Prof.  Larry  Crum 
and  co-workers.  In  addition  to  reaching  an  understanding  of  a  puzzling  “uni¬ 
versality”  of  rain  noise  at  low  rainfall  rates,  we  have  for  the  first  time  pro¬ 
duced  entirely  numerical  rain  noise  spectra  resembling  the  observed  ones. 
As  a  consequence  of  the  publicity  reached  by  this  work,  the  PI  was  asked 
to  provide  a  review  paper  on  the  problem  for  the  Annual  Reviews  of  Fluid 
Mechanics. 

Work  on  the  second  problem  has  been  based  on  the  initial  observation 
of  the  PI  (and,  independently  of  Dr.  W.  Carey),  originally  formulated  in  a 
paper  presented  in  1985  at  the  Nashville  meeting  of  the  Acoustical  Society  of 
America,  that  bubble  clouds  could  oscillate  in  collective  modes  at  frequencies 
as  low  as  a  few  tens  of  Hz,  even  though  the  constituent  bubbles,  in  isolation, 
might  have  natural  frequencies  of  tens  of  kHz  or  more.  This  hypothesis 
has  now  been  widely  confirmed  both  in  the  laboratory  and  in  the  field  and 
constitutes  the  accepted  answer  to  a  puzzle  that  had  remained  unsolved  since 
the  original  measurements  of  Knudsen  during  World  War  H. 

The  results  of  the  work  are  documented  in  a  series  of  papers  listed  below. 
The  first  page  of  each  paper  is  also  reproduced  below.  Additional  information 
on  the  work  carried  out  under  this  grant  is  provided  in  the  annual  summaries 
submitted  to  ONR  and  published  in  part  in  the  series  of  annual  volumes  titled 
Ocean  Acoustics  Program  Summary. 
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X' generalization  of  the  impulse  and  virial  theorems 
with  an  application  to  bubble  oscillations 

By  HASAN  N.  OGUZ  and  ANDREA  PROSPERETTI 

Department  of  Mechanical  Engineering  The  Johns  Hopkins  University. 

Baltimore,  MD  21218,  uSA 

(Received  11  August  1989  and  in  revised  form  29  January  1990) 

Tn  the  first  Dart  of  the  paper  it  is  shown  that  the  impulse  and  virial  theorems  of 

‘ "n  the  epptcxi.ete  study 

on  the  basis  of  the  linear  theory  of  Bjerknes  forces. 


1.  Introduction  .  ^  •  j  u 

following  relation :  ,  v  ^  r 

I  ^5ndSi=  [lu-«n-(n-u)«]d5B.  (El) 

SX  1986)  have  given  a  number  of  examples  ot  the  appheat.on  of  th.s  theorem 

“XtrtntegXheorem,  only  valid  tor  an  )»«"«' STvious 
Benjamin  (1987)  and.  in  a  different  way,  by  Longuet-Higgms  (1989).  In  the  previo 

notation  this  theorem  may  be  wTitten 

-Sp  f  =  f  (2>-7>J(^-”)d^i-  (^-2) 

”  Imthe‘present  paper  ue  shall  generalise  these  results  in  two  direetions.  First,  we 


Engineering  Science,  Fluid  Dynamics  (pp.  191  202) 

©  1990  World  Scientific  Publishing  Company 


BUBBLE  ENTRAPMENT  BY 
AXISYMMETRIC  CAPILLARY  WAVES 


HASAN  N.  06UZ  and  ANDREA  PROSPERETTI 

Department  of  Mechanical  Engineering 
The  John!  Hopkins  University 
Baltimore,  MD  ilZlS 


ABSTRACT 

Experiment  and  computation  show  that  liquid  drops  impacting  on  a  liq¬ 
uid  surface  can  entrain  air  bubbles  which  play  a  very  important  role  in 
the  underwater  noise  of  rain.  To  explore  the  physical  mechanism  respon¬ 
sible  for  this  process,  the  dynamics  of  a  large-amplitude,  circular  capillary 
wave  is  studied  numericaUy.  It  is  found  that,  as  the  wave  propagates  to¬ 
ward  the  axis  of  symmetry,  nonlinear  effects  can  lead  to  a  steepening  and 
overturning  strikingly  similar  to  those  found  in  the  impacting  drop  case. 
When  it  reaches  the  axis  of  symmetry,  in  suitable  conations,  the  wave  is 

also  able  to  entrap  air  bubbles.  •  '  . 


1  .Introduction 

The' fact  that  liquid  drops  falling  on  a  liquid  surface  may  entrain  air  bubbles  has 
been  known  for  some  time.*  However,  it  is  only  very  recently  that  the  detailed 
mechanics  of  this  process  has  been  elucidated  in  a  beautiful  set  of  experiments 
by  Pumphrey  and  Crum.*-='  Among  other  results,  these  authors  have  shown  that 
bubble  entrapment  occurs  in  a  small  and  well-defined  region  of  the  parameter 
space  formed  by  the  radius  R  of  the  impacting  drop  and  its  impact  velocity  V 
(Fig.  1).  Furthermore,  they  have  also  shown  that  very  little  acoustic  energy  is 
radiated  into  the  liquid  by  the  impact  itself,  the  most  significant  noise  source 
being  the  volume  pulsations  of  the  entrained  bubbles.  ^ 

These  observations  have  a  very  important  bearing  on  the  underwater  noise 
produced  by  rain.  The  steep  line  on  the  left  of  Fig.  1  shows  the  terminal  velocity 
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BUBBLE  ACTIVITY  INDUCED  BY  HIGH-POWER  MARINE  SOURCES 

G.B.  Cannelii’,  E.  D’Ottavi*,  and  A.  Prosperetti*' 

*CNR,  Istituto  di  Acustica  ”O.M.  Corbino”,  Via  Cassia  1216,  00189  Roma,  Italy. 
‘‘Department  of  Mechanical  Engineering,  The  Johns  Hopkins  University, 
Baltimore,  MD  21218,  U.S.A. 


Abstract:  This  paper  presents  some  preliminary  results  of 
a  systematic  study  of  the  bubble  phenomena  associated  with 
high-power  sparker-based  acoustic  marine  sources  of  a  novel 
type.  Data  taken  at  sea  and  in  a  laboratory  tank  are  pre¬ 
sented  and  an  attempt  is  made  to  interpret  them  in  the  light 
of  a  simple  theoretical  model. 


Introduction 

The  function  of  any  marine  source  of  seismic  energy  is  to 
generate  a  powerful  acoustic  wave  in  the  water  which  is  em¬ 
ployed  to  explore  sub-bottom  sea  structures.^  As  monopoles 
are  the  most  efidcient  acoustic  sources  in  a  fluid,  it  is  natural 
that  bubble  phenomena  play  a  prominent  role  in  the  func¬ 
tioning  of  most  of  the  sources  developed  to  date.  A  drawback 
is  however  the  fact  that,  in  most  cases,  the  bubbles  produced 
by  the  sources  give  rise  to  more  than  one  pressure  pulse. 
The  pressure  waves  associated  to  the  bubble  implosion  and, 
.  in  some  cases,  to  the  following  oscillations,  contaminate  the 
signal  and  require  sophisticated  data-processing  techniques 
for  their  removal.^  The  need  to  understand  these  bubble  pro¬ 
cesses  becomes  more  and  more  pressing  with  the  contempo¬ 
rary  trend  toward  the  use  of  more  powerful  acoustic  sources 
to  penetrate  deeper  into  the  sub-bottom  or  to  propagate 
underwater  signals  over  large  distances  for  long-range  com¬ 
munication. 

These  considerations  motivate  a  systematic  study  cur¬ 
rently  under  way  of  the  bubble  activity  induced  by  the  high- 
power,  paraboloidal,  sparker-based  marine  source  recently 
developed  at  the  Istituto  di  Acustica  ”O.M.  Corbino”.^’^  In 
comparison  with  traditional  devices,  this  new  source  is  char¬ 
acterized  by  a  more  controllable  electronics  and  mechanics 
which  makes  it  considerably  more  versatile.  A  set  of  experi¬ 
mental  tests  both  in  the  laboratory  and  at  sea  has  been  initi¬ 
ated  to  study  the  dependence  of  its  performance  on  different 
parameters  and  experimental  conditions.  In  this  paper  we 
present  some  preliminary  results  of  this  study,  in  particular 
concerning  the  effect  of  the  liquid  temperature.  On  the  ba- 
j  sis  of  these  experimental  results,  a  simple  theoretical  model 
'  is  developed  to  give  a  first  interpretation  of  the  observed 
f  phenomena. 


Experimental  tests 

Some  preliminary  tests  were  carried  out  at  sea  off  the 
coast  of  Fiumicino  (Rome)  in  order  to  compare  the  behav¬ 
ior  of  the  paraboloidal  sparker-based  source  with  that  of  a 
conventional  sparkarray.  The  principle  used  to  produce  the 
pressure  pulses  in  the  two  sources  is  the  same,  namely  the 
discharge  between  pairs  of  electrodes  in  the  liquid  of  a  large 
electrical  charge  stored  in  a  capacitor  bank.  However,  in  the 
paraboloidal  source,  the  electric  spark  gap  is  positioned  at 
the  focus  of  a  metallic  paraboloidal  reflector  and  the  pulse 
generated  is  characterized  by  very  high  power  in  a  wide  range 
of  low  and  middle  frequencies  (0.1  to  15  kHz)  not  usually 
available  in  standard  commercial  devices.  ' 


Figure  1:  Comparison  bettneen  the  acoustic  signatures  gen¬ 
erated  by  the  paraboloidal  sparker-based  source  (a)  and  of  a 
traditional  sparkarray  (b)  in  the  same  experimental  condi¬ 
tions  in  sea  water.  The  firing  energy  is  300  J. 

Figure  1  shows  a  comparison  between  the  signatures  of 
the  two.  sources  for  a  firing  energy  of  300  J.  [Note  that,  in 
this  as  in  the  following  figures,  the  polarity  is  reversed  so 
that  a  conipression  appears  as  a  decreeise  of  the  amplitude. 
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Bubble  entrainment  by  the  impact  of  drops  on 

liquid  surfaces 

By  HASAN  N.  OGUZ  and  ANDREA  PROSPERETTI 

Department  of  Mechanical  Engineering,  The  Johns  Hopkins  University, 

Baltimore,  MD  21218,  USA  . 

(Received  19  ^lay  1989  and  in  revised  form  2  March  1990) 

The  impact  of  a  drop  on  the  plane  surface  of  the  same  liquid  is  studied  numerically. 
The  accuracy  of  the  calculation  is  substantiated  by  its  good  agreement  with  available 
experimental  data.  An  attempt  is  made  to  explain  the  recent  observation  that,  in  a 
restricted  range  of  drop  radii  and  impact  velocities,  small  air  bubbles  remain 
entrained  in  the  liquid.  The  implications  of  this  process  for  the  underwater  sound  due 
to  rain  are  considered.  The  numerical  approach  consists  of  a  new  formulation  of  the 
boundary-element  method  which  is  explained  in  detail.  Techniques  to  stabilize  the 
calculation  in  the  presence  of  strong  surface-tension  effects  are  also  described. 


1.  Introduction 

The  fact  that  a  water  drop  impinging  on  a  water  surface  may  lead  to  the 
entrapment  of  an  air  bubble  at  the  bottom  of  the  crater  that  it  produces  has  been 
known  for  some  time  (Franz  1959).  However,  the  intricate  features  of  this 
phenomenon  have  only  recently  been  clarified  by  Pumphrey  &  Crum  (1988,  see  also 
Pumphrey,  Crum  &  Bjorno  1989)  in  the  course  of  an  experimental  study  of  the 
mechanism  of  rain  noise.  While  their  results  will  be  summarized  in  some  detail  in  the 
next  section,  we  mention  here  that  they  discovered  that,  far  from  being  a  random 
event  as  suggested  by  Franz,  bubbles  remain  entrapped  under  very  well-defined 
conditions  and  lead  to  a  very  substantial  noise  emission  in  the  water.  The  significant 
consequences  of  this  finding  for  the  generation  of  underwater  noise  by  rain  have  been 
addressed  elsewhere  (Prosperetti,  Pumphrey  &  Crum  1989)  and  are  summarized  in 
the  last  section  of  this  paper.  Here  we  wish  to  study  theoretically  the  fluid  dynamics 
of  crater  formation  and  bubble  entrapment  and  develop  a  qualitative  understanding 
of  the  physics  involved  in  these  processes.  Our  main  tool  is  numerical,  and  the 
dominant  effect  of  surface  tension  on  the  process  has  required  the  development  of  a 
stable  boundary-integral  method  which  is  of  interest  in  itself  and  is  described  in 

detail.  .  u  i  + 

The  study  of  drop  impact  on  liquid  surfaces  has  a  long  history  that  goes  back  to 

the  end  of  the  last  century  when  Worthington  (1894)  studied  the  process  by  means 
of  single-flash  photography.  Quite  famous  are  also  the  high-speed  photographs  of 
Edgerton  (Edgerton  &  Killian  1939).  These  researchers,  as  well  as  others  such  as 
Engel  (1967)  and  ilacklin  &  Metaxas  (1976),  studied,  however,  impacts  at  relatively 
high  speeds  and  therefore  missed  many  interesting  aspects  of  the  process  that  were 
recently  identified  by  Pumphrey  and  Crum.  One  of  the  early  applications  of  the  MAC 
code  also  dealt  with  drop  impact  on  a  liquid  layer  (Harlow  &  Shannon  1967 o,  6),  but 
again,  as  pointed  out  by  Carroll  &  Mesler  (1981  o,  b),  many  subtleties  were  missed 
because  of  the  incomplete  treatment  of  surface-tension  phenomena.  The  same 


Bubble  oscillations  in  the  vicinity  of  a  nearly  plane  free  surface 

Hasan  N.  Oguz  and  Andrea  Prosperetti 

Department  of  Mechanical  Engineering,  The  Johns  Hopkins  University,  Baltimore,  Maryland  21218 

(Received?  August  1989;  accepted  for  publication  7  December  1989) 

The  linear  oscillation  frequency  of  a  bubble  in  the  vicinity  of  a  distorted  plane  free  surface  is 
calculated  by  a  perturbation  method.  The  approximate  expression  found  is  compared  with 
numerical  results  valid  for  surface  deformations  of  arbitrary  magnitude.  It  is  found  that  the 
approximate  analytical  result  is  quite  good,  provided  that  the  deformation  is  small  compared 
with  the  depth  of  immersion  of  the  bubble.  It  is  also  shown  that,  unless  the  deformation  of  the 
free  surface  extends  to  distances  at  least  of  the  order  of  an  acoustic  wavelength,  the  “image’’ 
bubble  has  the  same  source  strength  of  the  real  bubble  so  that  a  dipolar  acoustic  emission  can 
be  expected  in  spite  of  the  deformation  of  the  surface. 

PACS  numbers:  43.35.Pt,  43.30.Lz,  43.30.Nb 


INTRODUCTION 

Recent  work  shows  that  a  substantial  contribution  to 
oceanic  ambient  noise  is  caused  by  the  free  oscillations  of  air 
bubbles  entrained  by  breaking  waves,  by  impacting  rain 
drops  or  sprays,  or  possibly  by  capillary  waves  of  limiting 
form.*"’  In  all  these  cases,  at  the  instant  when  the  surface 
closes,  the  resulting  bubble  is  not,  in  general,  in  equilibrium. 
The  excess  initial  energy  is  dissipated  in  the  course  of  shape 
and  volume  oscillations  that  give  rise  to  the  acoustic  emis¬ 
sion. 

Although  a  great  deal  is  known  about  the  oscillations  of 
bubbles  in  unbounded  liquids  (see,  e.g..  Refs.  8-13),  it  is  not 
obvious  that  this  information  is  applicable  to  the  circum¬ 
stances  described  above  because,  during  the  acoustic  emis¬ 
sion,  the  bubble  is  near  the  ocean  surface,  and  this  surface  is 
not  plane, 

A  typical  example  is  shown  in  Fig.  1  (with  the  permis¬ 
sion  of  H.  C.  Pumphrey  and  L.  A.  Crum),  which  is  a  frame 
from  a  high-speed  movie  film  of  a  water  drop  falling  on  an 
undisturbed  water  surface.  The  left  side  of  the  frame  shows  a 
bubble  that  has  just  detached  from  the  bottom  of  the  crater 
created  by  the  impact.  The  right  side  of  the  frame  shows  the 
oscilloscope  trace  produced  by  a  hydrophone  driven  by  the 
sound  pulse  emitted  by  the  bubble.  The  frequency  of  the 
damped  sinusoid  measured  from  this  photograph  is  approxi¬ 
mately  7  kHz,  and  the  bubble  diameter,  measured  from  a 
later  frame  in  which  the  bubble  is  in  equilibrium,  is  0,95  mm. 
It  is  most  remarkable  that,  in  spite  of  the  proximity  of  the 
highly  distorted  free  surface,  this  frequency  is  in  almost  per¬ 
fect  agreement  with  the  natural  frequency  of  a  bubble  of  the 
same  radius  in  an  unbounded  fluid,  which  equals  6.78  kHz. 

In  the  present  paper,  we  try  to  find  a  theoretical  justifi¬ 
cation  for  this  fact  on  the  basis  of  a  perturbative  solution  of 
the  problem.  The  approximation  results  are  checked  against 
numerical  ones,  and  their  domain  of  validity  is  established. 
In  the  last  section,  we  offer  some  comments  on  the  expected 
acoustic  radiation  from  the  bubble  under  these  conditions. 

The  basis  of  our  approximation  consists  in  the  fact  that, 
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in  many  situations  of  interest,  the  distortion  of  the  free  sur¬ 
face  has  a  scale  larger  than  the  bubble  radius  and  the  distance 
of  the  bubble  from  the  free  surface  is  not  small.  The  first 
circumstance  arises  because,  in  general,  the  events  that  ren¬ 
der  the  free  surface  multiply  connected  are  so  energetic  as  to 
produce  disturbances  on  a  much  larger  scale  than  the  bub¬ 
ble.  Second,  in  the  neighborhood  of  the  point  from  which  the 
bubble  detaches,  the  curvature  is  very  large.  As  soon  as  the 
bubble  is  formed,  therefore,  under  the  action  of  surface  ten¬ 
sion,  this  local,  high-curvature  deformation  disappears  very 


FIG.  1.  A  frame  from  a  high-speed  cinematographic  sequence  showing  the 
entrapment  of  an  air  bubble  when  a  drop  hits  a  liquid  surface  (courtesy  of 
H.  C.  Pumphrey  and  L.  A.  Crum).  The  left  half  of  the  frame  shows  the 
physical  event.  One  can  see  the  free  surface  indented  by  the  crater  caused  by 
the  impact  and,  directly  under  it,  the  bubble  that  has  just  detached.  The 
shadow  in  the  lower  right  is  a  hydrophone  that  drives  an  oscilloscope  the 
trace  of  which  is  shown,  vertically,  in  the  right  half  of  the  frame.  The  right 
half  of  the  picture  shows  therefore,  in  a  sense,  the  “soundtrack”  of  the  left 
half.  The  damped  sinusoid  of  the  oscilloscope  trace  is  the  acoustic  signal 
produced  by  the  volume  oscillations  of  the  bubble.  A  study  of  the  angular 
distribution  of  the  signal  reveals  the  expected  dipole  pattern.  See  Refs.  3-5 
for  further  details. 
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Abstract— 'Ry  means  of  an  effeclive  equation  model  for  the  propaga¬ 
tion  of  pressure  waves  in  a  bubbly  liquid,  the  normal  modes  of  oscilla¬ 
tion  of  regions  of  bubbly  liquid  in  an  otherwise  pure  liquid  are  calcu¬ 
lated  for  some  simple  geometries.  It  is  shown  that  the  frequencies  of 
oscillation  of  such  bubble  clouds  can  be  very  much  lower  than  those  of 
the  constituent  bubbles  in  isolation  and  fall  well  within  the  range  where 
substantial  wind-dependent  noise  is  observed  in  the  ocean.  A  compari¬ 
son  with  some  experimental  data  very  strongly  supports  the  theoretical 
results. 

I.  Introduction 

IN  THE  OCEAN,  bubble  clouds  are  incessantly  formed  in 
the  upper  layers  by  the  breaking  of  waves  and  are  trans¬ 
ported  to  depths  of  tens  of  meters  by  Langmuir  circulation, 
turbulence,  and  other  mechanisms  [1],  [2].  One  might  expect 
these  bubbles  to  contribute  substantially  to  the  underwater 
ambient  noise,  since  wave  breaking  is  a  catastrophic  event 
that  will  certainly  impart  them  an  appreciable  initial  energy 
[3]_[7].  That  a  substantial  part  of  this  acoustic  emission 
could  be  at  frequencies  well  below  a  few  kilohertz,  however, 
may  at  first  sight  be  surprising.  Indeed,  from  the  approximate 
relation 


relating  the  bubble  radius  a  to  its  natural  frequency  Vq  = 
aijlit,  one  finds,  for  example,  Vq  =  2.8  kHz  for  a  bubble 
having  the  radius  of  1  mm  in  water  (density  p  =  1  g/cm^)  at 
normal  ambient  pressure  {Pq  -  1  atrn).  However,  bubbles  in 
a  cloud  are  in  effect  coupled  oscillators,  and  one  can  there¬ 
fore  expect  the  existence  of  normal  modes  of  the  cloud  itself 
at  a  substantially  lower  frequency  than  that  of  the  individual 
constituent  bubbles.  An  alternative  argument  leading  to  the 
same  conclusion  may  be  based  on  the  observation  that  the 
speed  of  sound  in  a  bubbly  mixture  can  be  an  order  of 
magnitude  smaller  than  in  the  pure  liquid,  even  at  gas-volume 
fractions  less  than  1%.  Therefore,  to  a  first  approximation,  a 
region  of  bubbly  liquid  can  be  considered  as  enclosed  by  a 
rigid  boundary  and  will  then  possess  normal  modes.  How¬ 
ever,  the  fact  that  the  boundary  is  not  really  rigid  has  the 
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consequence  that  the  acoustic  energy  “trapped”  in  the  bub¬ 
bly  region  will  leak  out  into  the  pure  liquid. 

The  frequencies  o}/.  of  these  normal  modes  are  readily 
estimated  as  follows  [7].  Let  the  cloud  have  linear  dimen¬ 
sions  of  order  L.  Then,  from  the  analogy  with  similar 
systems  [8],  one  expects  eigenfrequencies  of  the  order  of 
oj*  =  Arc„/2  L ,  A:  =  1 , 2 ,  •  •  • ,  where  c„  is  the  effective  speed 
of  sound  in  the  bubbly  mixture  given  approximately  by  the 
well-known  expression  [9]-[ll]: 


where  (3  is  the  gas-volume  fraction.  If  the  constituent  bubbles 
have  a  radius  of  order  d,  the  gas-volume  fraction  may  be 
estimated  by  |3  ~  {ayL^)N,  where  N  is  the  number  of 
bubbles  in  the  cloud.  In  this  way  we  find: 


(jif.  l^ 


(3) 


The  volume  fraction  here  is  raised  to  such  a  low  power  that 
the  result  is  nearly  independent  of  this  variable.  Hence  all  the 
frequencies,  and  in  particular  the  lowest  one  corresponding  to 
k  =  1 ,  are  predicted  to  decrease  essentially  as  the  cube  root 
of  the  bubble  number.  For  N  ~  1000,  the  frequency  reduc¬ 
tion  with  respect  to  the  single-bubble  case  would  be  of  one 
order  of  magnitude  [7]. 

The  preceding  estimate  can  be  put  on  a  firmer  ground  by 
means  of  a  model  of  a  bubbly  liquid  which  is  adequate  up  to 
gas-volume  fractions  of  a  few  percent.  The  model,  which 
regards  the  gas-liquid  mixture  as  a  continuum  governed  by 
effective  equations,  is  essentially  due  to  Foldy  [12],  and  has 
been  rederived  by  many  others,  in  particular  in  [13].  Its- 
predictions  have  been  compared  with  experiment,  most  re¬ 
cently  in  [14],  and  a  very  good  agreement  was  found.  For 
completeness  we  briefly  present  the  linearized  version  of  this 
model  in  the  next  section.  A  series  of  papers  containing 
detailed  applications  of  it  to  the  generation  and  scattering  of 
sound  by  bubble  clouds  of  different  shape  and  in  a  number  of 
situations  is  being  prepared.  The  present  study  is  a  survey  of 
the  results  on  the  generation  of  noise  contained  in  these 
papers. 


II.  Mathematical  Formulation 

The  averaged  continuity  and  momentum  equations  applica¬ 
ble  to  a  bubbly  liquid  at  small  gas-volume  fraction  may  be 
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A  number  of  processes  in  which  air  is  entrained  in  a  flow  appear  to  involve  the 
formation  of  a  thin  air  film  between  a  relatively  fast  liquid  stream  and  a  region  of 
slow  recirculation.  Eventually,  the  film  breaks  into  bubbles.  This  study  addresses  a 
possible  mechanism  causing  this  process.  The  linear  stability  of  a  vertical  film  of  a 
viscous  gas  bounded  by  liquid  in  uniform  motion  on  one  side,  and  by  liquid  at  rest 
on  the  other  side,  is  studied.  Instabilities  are  found  that,  depending  on  the  parameter 
values  of  the  undisturbed  flow,  are  controlled  by  two  basic  mechanisms.  One  is  due 
to  the  velocity  jump  across  the  film  and  can  be  related  to  the  usual  Kelvin-Helmholtz 
instability.  The  second  one  is  controlled  by  the  viscosity  jump  across  the  air-liquid 
interfaces.  The  relation  between  the  remainder  of  the  discrete  spectrum  and  the 
spectrum  of  other  parallel  shear  flows  bounded  by  solid  or  free  surfaces  is  also 
discussed. 


1.  Introduction 

The  entrainment  of  air  in  a  flow  is  an  important  process  very  frequently 
encountered.  The  ecological  balance  of  water  bodies,  from  small  lakes  to  entire 
oceans,  is  critically  dependent  on  the  amount  of  dissolved  o.xygen.  Aeration  is  a 
standard  technique  of  water  treatment.  Furthermore,  the  formation  and  detachment 
of  bubbles  is  an  inherently  noisy  process  to  which  much  of  the  oceanic  ambient  noise 
over  a  large  frequency  range  from  hundreds  of  Hz  to  many  tens  of  kHz  can  be 
ascribed.  In  spite  of  this  widespread  occurrence,  not  much  seems  to  be  known  about 
the  basic  mechanisms  by  which  entrainment  takes  place.  In  a  paper  devoted  to  air 
entrainment  in  a  wave  breaking  in  the  spilling  mode,  Longuet-Higgins  &  Turner 
(1974)  mention  the  ‘over-running  of  air  by  the  advancing  front’  of  water  and  the 
‘self-aeration’  of  thin,  highly  turbulent  flows  which  develops  when  the  turbulent 
boundary  layer  on  the  bottom  reaches  the  surface.  While  certainly  correct  and 
adequate  for  the  purposes  of  their  study,  these  statements  are  rather  vague  as  to  the 
precise  nature  of  the  process.  A  literature  search  has  not  produced  much  more 
detailed  information  than  this.  In  the  present  paper  we  wish  to  investigate 
theoretically  a  possible  mechanism  by  which  air  can  be  entrained  in  flows.  Although 
not  the  only  one,  this  mechanism  appears  to  be  of  sufficiently  widespread  occurrence 
to  warrant  its  investigation. 

A  consideration  of  several  examples  of  entraining  flows  suggests  that  a  possible 
mechanism  involves  the  development  and  instability  of  a  thin  air  film  at  the 
boundary  between  two  liquid  currents.  The  clearest  example  of  this  process  is  offered 
by  a  jet  falling  into  a  liquid  pool.  In  a  high-viscosity  liquid  an  air  film  surrounding 
the  jet  can  be  clearly  discerned  for  several  diameters  below  the  free  surface  (Lin  & 
Donnelly  1966).  The  film  develops  a  wavy  structure,  the  amplitude  of  which 
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When  raindrops  with  a  diameter  of  the  order  of  1  mm  hit  a  ^ 

entrain  air  bubbles  that  radiate  noise  in  the  course  of  volume  oscillations.  The  paper 
presents  a  model  of  the  underwater  noise  of  rain  produced  by  this  process.  The  dep 
of  submergence,  radius,  and  initial  energy  of  the  entrained  bubbles  are  obtamed 
numerically  for  a  number  of  drop  sizes.  The  bubbles  are  assumed  to  radiate  as 
Ses  and  the  total  underwater  noise  is  calculated  by  intepating  over  the  size  of 
the  entwining  rain  drops.  The  results  are  compared  both  vv^.th 

of  single-drop  impacts  and  field  data  of  ram  noise.  It  is  found  that  the  model  gives 
somewhat  larger  bubbles  than  are  observed  experimentally.  As  a  consequ^ence,  e 
characteristic  spectral  peak  of  rain  is  predicted  to  occur  at  a  somewhat  lo\ 
frequency  than  found  in  experiment.  However  the  level  of  the  Pf^ 
agreement  with  data.  The  amount  of  noise  due  to  the  process 

is  also  estimated  and  found  to  be  several  orders  of  magnitude  lov  er  than  the  data. 
Therefore,  in  spite  of  some  deficiencies  of  the  model  and  of  the  computational  results, 
the  proposed  mechanism  for  the  underwater  noise  of  ram  is  strongly  supported  by 

this  study. 


1.  Introduction 

The  underwater  noise  produced  by  rain  has  recently  been  measured  by  a  number 
of  groups  and  in  different  conditions  (Scrimger  198o;  Scrimger  et  al.  ^  ’ 

Evins  &  Yee  1989;  Nystuen  1986;  Nystuen  &  Farmer  1989;  Pumphrey,  Crum  <L 
Bjorno  1989).  A  striking  finding  common  to  all  these  studies  is  e  presen  y 

prominent  and  well-defined  spectral  peak  at  a  frequency  of  e 

Quite  unexpectedly,  the  position  and  general  shape  of  this  peak  are  found  to 
Lependen?  of  the  rainfall  rate,  and  even  of  the  size  djlnbution  of  the  ram  drops^ 
A  striking  example,  from  Scrimger  et  al.  (1987),  is  repoduced  in  figure  1.  Fig«r^) 
shows  the  measured  raindrop  size  distributions  for  two  ram  eventj 
noise  of  which  is  shown  in  figure  1  (6).  It  can  be  seen  that  the  size 
,  .very  different,  while  the  measured  spectra  are  strikingly  alike.  Some  other  examples 

can  be  found  in  Scrimger  ef  ah  (1987).  an 

In  a  recent  study  (Prosperetti,  Pumphrey  &  Crum  1989)  we  have  Foposed  an 
explanation  for  these  findings  based  on  the  fluid  mechanics  of  the  impact  of  ^P 
liquid  surfaces.  The  starting  point  is  work  by  Pumphrey  &  Crum  ^"“P  .  ^ 

aZ.  (1989),  Pumphrey  &  Crum  (1990)  and  Pumphrey  &  Elmore  (1990)  ^cb 
presents  conclusive  evidence  that  most  of  the  noise  produce  y  impac  i  g  . 

at  normal  incidence  is  not  actually  produced  by  the  impact  itself,  b"t  bj  ^^aU  a  " 
bubbles  entrained  by  the  impacting  drop.  A  remarkable  feature  of  this  process 
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. ,  It  is  well  known  that  ocean  ambient  noise  levels  in  the  frequency  range  from  a  few  hundred 
hertz  to  several  tens  of  kilohertz  are  well  correlated  with  wind  speed.  A  physical  mechanism 
-  that  could  account  for  some  of  this  sound  generation  is  the  production  of  bubble  clouds  by 
breaking  waves.  A  simple  laboratory  study  of  the  sound  generated  by  a  column  of  bubbles  is 
reported  here.  From  measurements  of  the  various  characteristics  of  this  column,  good  evidence 
is  obtained  that  the  bubbles  within  the  column  are  vibrating  in  a  collective  mode  of  oscillation. 

Based  upon  an  assumption  of  collective  oscillations,  analytical  calculations  of  the  predicted 
frequency  of  vibration  of  this  column  as  well  as  the  dependence  of  this  frequency  on  such 
parameters  as  bubble  population  and  column  geometry  agree  closely  with  the  measured  values. 

These  results  give  evidence  that  the  bubble  plumes  generated  by  breaking  waves  can  be  a 
strong  source  of  relatively  low  frequency  (  <  1  kHz)  ambient  noise. 

PACS  numbers:  43.30.Nb  . 


INTRODUCTION 

The  literature  contains  substantial  evidence  that 
oceanic  ambient  noise  is  wind-dependent  down  at  least  to 
frequencies  of  a  few  hundred  hertz.’”*^  The  nature  of  the 
source  of  these  low-frequency  acoustic  emissions  is,  how¬ 
ever,  still  obscure.  In  view  of  this  seeming  correlation  with 
breaking  waves,  it  was  suggested  independently  by  Carey**'^ 
and  ProsperettP"*^  that  the  collective  oscillations  of  bubble 
clouds  could  explain  these  emissions.  It  is  the  purpose  of  the 
present  paper  to  present  an  experimental  confirmation  of 
this  possibility,  and  to  show  that  analytical  predictions  agree 
closely  with  the  measurements. 

It  is  well  known  that  bubbles  abound  in  the  first  several 
meters  under  the  ocean  surface.  Their  origins  are  varied  (im¬ 
pact  of  splashes,  sprays,  rain  drops,  biological  activity,  capil¬ 
lary  waves  of  limiting  form,  and  other  mechanisms),  but 
probably  the  most  significant  source  of  these  bubbles  is  the 
breaking  of  surface  waves. The  precise  nature  of  the 
mechanism  of  air  entrainment  in  a  breaking  wave  is  still 
poorly  understood,  but  it  is  plausible  that  at  the  moment  at 
which  the  liquid  surface  closes  on  itself  to  entrap  bubbles,  its 
local  velocity  is  nonzero.  The  bubble  created  in  this  way  has 
therefore  some  initial  kinetic  energy  and,  unless  it  is  spheri¬ 
cal  and  with  a  balanced  initial  pressure,  some  initial  poten¬ 
tial  energy  as  well.  Since  a  bubble  may  roughly  be  regarded 
as  a  mechanical  oscillator,  one  expects  this  initial  energy  to 
give  rise  to  volume  pulsations  with  associated  acoustic  emis¬ 
sions.  The  reason  why  such  a  simple  picture  is  unacceptable 
as  an  explanation  for  the  observed  low-frequency  noise  is 
that  the  natural  angular  frequency  of  a  bubble  of  radius  a 
is  given  approximately  by 

00  =  {\/a){lypa/py'^, 

“’On  leave  from  the  Sung  Kyun  Kwan  University,  Seoul,  Republic 
of  Korea. 


where /?o  is  the  ambient  pressure,  p  is  the  water  density,  and 
y  is  the  ratio  of  specific  heats.  With  /?(,  —  lO'*^  Pa,  p  —  10*^ 
kg/m’\  y—  1.4,  this  relation  would  require  a  radius  of  about 
7  mm  for  a  frequency / ^  =  500  Hz,  and  it  is  by  no 

means  clear — and  actually  rather  unlikely — that  such  large 
bubbles  are  created  in  significant  numbers. 

The  collective  oscillation  hypothesis  circumvents  this 
problem  by  appealing  to  the  well-known  fact  that  a  system  of 
coupled  oscillators  (such  as  the  bubbles  in  the  cloud,  for 
which  the  coupling  is  provided  by  mutual  hydrodynamic 
and  acoustic  interaction)  possesses  normal  modes  the  fre¬ 
quencies  of  which  can  be  substantially  lower  than  the  natural 
frequencies  of  the  individual  oscillators.  An  intuitive  argu¬ 
ment  leading  to  the  same  conclusion  is  as  follows.  It  is  well 
known  that  a  mixture  of  air  bubbles  and  liquid  has  a  sound 
speed  much  lower  than  that  of  the  pure  liquid  even  at  gas 
volume  fractions  as  low  as  1%.  As  a  rough  approximation, 
one  can  therefore  consider  a  bubble  cloud  as  an  acoustic 
medium  surrounded  by  a  rigid  enclosure,  a  system  which  is 
capable  of  normal  modes  of  oscillation  depending  on  its  lin¬ 
ear  dimensions.  Since  the  surrounding  liquid  is  not,  how¬ 
ever,  a  rigid  enclosure,  the  energy  trapped  in  the  cloud  will  ^ 
“leak  out”  and  be  detectable  as  acoustic  waves  away  from 
the  bubble  cloud.  A  very  simple  argument  shows  that  the 
ratio  of  the  minimum  cloud  eigenfrequency  (o„^  to  the  natu¬ 
ral  frequency  co^  of  the  constituent  bubbles  (assumed  to  be 
equal)  is  of  the  order  of 

where  [3  is  the  gas  volume  or  void  fraction  in  the  cloud  and 
Niy  is  the  total  number  of  bubbles.  For  example,  a  region  with 
linear  dimensions  of  the  order  of  10  cm  and  a  void  fraction  of 
1%  could  contain  enough  1  mm  radius  bubbles  to  cause  a 
frequency  reduction  of  one  order  of  magnitude. 

The  theoretical  literature  on  bubble  cloud  collective-os¬ 
cillation  phenomena  has  seen  several  contributions  in  the 
last  few  years.'^^’^'^*^  However,  so  far,  no  direct  experimen- 
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Abstract  Historical  Background 

Precipitation  plays  an  important  role  in  the  climate  of  a 
particular  region  and.  more  generally,  in  our  total  global 
environmenL  Its  occurrence  has  relevance  to  latent  heat  re¬ 
lease,  upward  mass  flux,  and  the  spatial  organization  of  con¬ 
vection.  Thus,  knowledge  of  the  rate  of  rainfall  and  its 
geographical  distribution  is  of  major  interest. 

Due  to  rain’s  spatial  and  temporal  discontinuity,  accurate 
measurements  are  difficult  even  over  land  because  rain  gauges 
at  discrete  sampling  sites  can  give  misleading  values.  Radar 
can  provide  a  more  complete  coverage,  but  its  general  use  is 
restricted  to  the  develop^  countries.  Rainfall  measurements 
over  the  ocean  where  -  according  to  estimates^  ~  80%  of  the 
earth’s  precipitation  occurs,  are  yet  more  difficult  to  obtain. 
Typically,  the  few  weather  stations  located  on  the  ocean  are 
on  islands,  a  circumstance  that  leads  to  an  inadequate  sampling 
of  the  global  distribution.  Shipboard  measurements  suffer 
problems  such  as  contamination  by  sea  spray,  platform  in¬ 
stabilities,  and  spatially-biased  sampling. 

Satellite-based,  remote  sensing  of  the  earth’s  environment 
probably  provides  the  highest  promise  for  adequate  global 
rainfall  monitoring.  Yet,  such  measurements  are  time-consum¬ 
ing,  costly,  and  suffer  the  important  limitation  of  “ground- 
truth,”  i.e.,  adequate  land-  or  sea-based  calibration. 

It  is  against  this  background  that  Walter  Munk’s  idea  of 
using  the  underwater  noise  produced  by  rain  to  monitor  rain¬ 
fall  rates  remotely  over  the  ocean  can  be  appreciated.  His 
suggestion  was  taken  up  by  Nystuen,  Scrimger,  and  others^'^ 
in  the  early  eighties  with  very  unexpected  and  intriguing  results. 


2  Nayal  Research  Reviews 


The  Office  of  Naval  Research  has  supported  a  research 
program  at  our  respective  institutions  for  the  past  few  years  in 
the  general  area  of  bubble-related  ambient  noise  in  the  ocean. 
In  the  course  of  our  study  of  this  topic,  we  came  to  discover 
some  interesting  things  about  the  underwater  sounds  produced 
by  the  impact  of  rain,  hail  and  snow  on  the  surface  of  the  sea. 
This  paper  describes  the  physical  processes  that  produce  these 
acoustic  emissions  and  shows  that  a  major  portion  of  precipi¬ 
tation  noise  is  caused  by  the  injection  of  gas  bubbles  into  the 
liquid  by  the  impact  of  rain  drops,  hail  stones  and  snow  flakes. 


Introduction 

The  distinctive  sound  of  a  raindrop  striking  a  puddle  of 
water  is  familiar  to  us  all.  It  is  one  of  those  unique  sounds  in 
nature  that  is  immediately  recognizable.  The  sound  of  a  pebble 
tossed  into  a  quiet  pond  has  a  similar  uniqueness  and,  to  most 
of  us,  even  evokes  a  sense  of  tranquillity.  The  soft,  acoustic 
environment  of  falling  snow  is  the  very  essence  of  stillness. 
We  have  examined  these  forms  of  precipitation  and  have  found 
that  the  noise  produced  when  they  impact  a  water  surface  has 
at  its  origin  an  oscillating  gM  bubble.  The  study  of  the  under¬ 
water  sounds  produced  by  precipitation  is  part  of  our  general 
study  of  bubble-related  ambient  noise  in  the  ocean,  which  has 
been  supported  by  the  Ocean  Acoustics  Program  of  the  Ocean 
Sciences  Directorate,  Office  of  Naval  Research. 
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ABSTRACT.  Experiment  and  theory  give  very  strong  indications  that  most  of  the  “"^erwa  e 
noise  of  rain  in  the  region  between  8-10  kHz  and  20-25  kHz  .s  due  to  the  radiations  of  bubbles 
entrained  by  the  rain  drops  impacting  on  the  ocean’s  surface.  We  present  a 

of  the  noise  spectrum  of  this  process  and  examine  the  effects  of  v.scos.ty  and  other  ^is 

result.  It  is  found  that  neglect  of  dissipative  effects  can  lead  to  bubble  sizes  greater  than  observed 

experimentally. 

1  Introduction 

Recent  experimental  work  by  Crum  and  co-workers  has  given  conclusive  ^ 

bubbles  entrained  when  a  liquid  drop  hits  a  plane  surface  are  acoustic  sources  of  greater 

magnitude  than  the  impact  itself  [1-3].  This  work  also  shows  that  for 

at  their  terminal  velocity,  only  a  restricted  range  of  drop  radii,  between  0.41  and  0.53 

mm  appro.ximately,  will  result  in  the  entrainment  of  a  bubble.  These  results  suggest  a 

possible  mechanism  for  the  peculiar  frequency  dependence 

that  seems  always  to  exhibit  a  spectral  peak  near  14  kHz  [  ,  ].  le  aroU 

the  acoustic  signal  depends  only  on  such  a  narrow  size  range  of  drop  ^ 

distributions  of  different  rain  events  would  be  very  similar  in  this  range, 

being  mostly  confined  to  the  total  drop  count.  Hence  the  acoustic  spectra  of  different  ram 

events  would  be  very  similar  in  relative  spectral  composition,  the  mam  difference  being 

the  absolute  “loudness”  of  the  noise  [6].  ,  n  , 

In  a  recent  study  [7]  we  have  attempted  to  reproduce  purely  theoretica  1^  without 
recourse  to  experimental  information,  the  underwater  noise  of  ram.  We 
and  the  initial  energy  of  the  bubbles  entrained  by  rain  drops,  and  assume  that  they  execute 
tod  hLonic  osciUations.  The  noise  spectra  generated  in  this  way  bear  a  striking 
resemblance  to  those  of  naturally  occurring  rain  noise,  except  that  the  peak 
occurs  at  a  frequency  of  8-9  kHz  rather  than  14  kHz.  In  the  Mowing  section  we  present 
a  summary  of  this  earlier  work,  and  then  devote  the  rest  of  the  paper  to  analysing  some 
of  the  possible  causes  of  this  discrepancy. 
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The  emission  and  scattering  of  sound  from  bubble  clouds  is  studied  thcoredcally.  It  is  shown 
that  clouds  having  a  size  and  air  content  similar  to  what  might  be  expected  as  a  consequence  of 
the  breaidng  of  ocean  waves  oscillate  at  frequencies  as  low  as  100  Hz  and  below.  Thus  cloud 
oscillations  may  famish  an  explanation  of  ’  the  substantial  amount  of  low-frequency 
wind-dependent  oceanic  ambient  noise  observed  experimentally.  Detailed  results  for  the 
backscattering  from  bubble  clouds — particularly  at  low  grazing  angles — are  also  presented  and 
shown  to  be  largely  compatible  with  oceanic  data  Although  the  cloud  model  used  here  is 
idealized  (a  uniform  hemispherical  cloud  under  a  plane  water  free-suiface),  it  is  shown  that  the 
results  are  relatively  robust  in  terms  of  bubble  size,  distribution,  and  total  air  content,  A  similar 
insensitiYiry  to  cloud  shape  is  found  in  a  companion  paper  [Sarkar  and  Prospererti,  J.  Acoust, 
Soc,  Am.  93,  3128-3138  0993)].  . 

"  PACS  numbers:  43.30.Nb,  43,30.Ft,  43.30.Lz 


INTRODUCTION 

Recent  research  has  brought  to  the  fore  the  impor¬ 
tance  of  bubble  clouds  at  the  ocean's  suriace  for  the  prop¬ 
agation  and  generation  of  underwater  sound  (Carey  and 
Bradley,  1985;  Carey  and  Browning,  1988;  Prosperetti, 
1985,  1988a,b;  Lu  et  at,  1990;  Yoon  e:  aL,  1991;  Lu  and 
Prosperetti,  1993).  On  the  one  hand,  it  has  been  realized 
that  the  clouds  can  oscillate  in  collective  modes  and  give 
rise  to  acoustic  emissions  as  low  as  a  few  tens  of  Hz  irre¬ 
spective  of  the  size  of  the  constituent  bubbles.  It  is  possible 
that  in  this  fact  an  e.xplanation  may  be  found. of  the  low- 
frequency  component  of  wind-dependent  oceanic  ambient 
noise  that  has  long^  puzzled  investigators  (Wenz,  1962; 
Kerman,  1984;  Carey  and  Wagstaff,  1986;  Kewiey  et  aL, 
1990).  Secondly,,  several  analyses  have  shown  that  the 
backscattering  produced  by  bubble  clouds  can  be  quite  sub¬ 
stantial  and  preliminary  estimates  mdicate  that  these  enti- 
ti«  may  . be  responsible  for  the  unexpectedly .  high  back- 
scattering  strengths  found  experim^^y. 

The  presence  of  bubbles  in  substantial  numbers  in  the 
uppermost  several  meters  of  the  ocean  surface  is  so  well 
known  that  a  few  references  will  be  sufficient  in  this  respect 
(Monahan,  1971;  Thorpe,  1982;  Farmer  and  Vagle,  1989). 
One  may- distinguish  between  the  regions  of  “fresh”  bub¬ 
bles,  newly  formed  by  breaking  waves,  and  the  “old”  bub¬ 
bles,  that  are  stabilized  by  still  unclear  mechanisms  and 
survive  long  after  their  formation.  Typically,  “fresh”  bub¬ 
bles  arc  relatively  highly,  concentrated  in  clouds  very  close 


‘^Present  address:  Applied  Physics  Labomtory,  Submarine  Technology 
Department,  Hydrodynamics  Group,  Johns  Hopkins  Road,  Laurel,  MD 
20723-^9, 
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jeon  305r^.  Korea. 
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to  the  sunace,  while  “old”  bubbles  are  much  less  dense  and 
are  transponed  downward  by  turbulence  and  Langmuir 
circulation.  We  refer  to  these  latter  agglomeration  as  bub¬ 
ble  “plumes”  reser/ing  the  term  “cloud”  for  the  former. 

In  this  paper  we  consider  both  the  active  and  passive 
facets  of  the  acoustic  behavior  of  bubble  clouds.  Our  model 
is  geometrically  simple  and  consists  of  hemispherical 
clouds  at  the  surface  of  a  plane  ocean.  However,  we  use  a 
relatively  complete  model  for  the  bubbly  region  that  goes 
beyond  previous  calculations  and  we  solve  exactly — rather 
than  in  an  approximate  fashion — the  emission  and  scatter¬ 
ing  problems.  In  our  basic  model  the  bubble  cloud  has  a 
uniform  spatial  distribution  of  gas  and  bubble  radiL  Later, 
we  consider  generalizations  of  this  model  and  try  to  derive 
conclusions  of  more  general  validity.  The  dependence  of 
the  results  on  the  cloud  shape  is  dearly  also  important,  and 
is  addressed  in  a  separate  publication  (Sarkar  and  Pros¬ 
peretti,  1993a). 

Crowthcr  (1980),  McDaniel  and  Gorman  (1982, 
1983),  and  McDaniel  ( 1987)  liave  applied  incoherent  scat¬ 
tering  methods  to  the  calculation  of  backscattering  from 
bubble  layers.  More  recently,  McDonald  •  (1991)  and 
Henyey  (1991)  have  treated  the  backscattering  from  bub¬ 
ble  plumes  by  the  Bom  approximation.  The  methods  used 
by  these  auffiors  are  not  suitable  for  the  denser  bubble 
assemblies  studied  in  this  paper.  A  detailed  comparison  of 
these  approaches  be  found  in  Sarkar  and  Prosperetti 
(1993a,  1993b). 

Perhaps  the  major  piece  of  information  missing  for  a 
quantitative  description  of  bubble  clouds  is  the  gas  concen¬ 
tration  by  Yoltime,  also  called  void  fraction  in  the  multi¬ 
phase  flow  literature.  Some  estimates  of  this  quantity  dur- 
mg  the  active  breaking  process  are  as  high  as  30% 
(Longuet-Higgins  and  Turner,  1974;  MeiviUc  et  aL,  1992), 
but  this  situation  of  extremely  large  void  fraction  must  be 
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INTRODUCTION 


There  will  be  but  few  of  my  readers  who  have  not,  in  some  heavy  shower  of  rain, 
beguiled  the  tedium  of  enforced  waiting  by  watching,  perhaps  half-unconsciously,  the 
thousand  little  crystal  fountains  that  start  from  the  surface  of  pool  or  river:  noting 
now  and  then  a  surrounding  coronet  of  lesser  jets,  or  here  and  there  a  bubble  that 
floats  for  a  moment  and  then  vanishes. 

It  is  to  this  apparently  insignificant  transaction,  which  always  has  been  and  always 
will  be  so  familiar,  and  to  others  of  a  like  nature,  that  1  desire  to  call  the  attention  of 
those  who  are  interested  in  natural  phenomena;  hoping  to  share  with  them  some  of 
the  delight  that  I  have  myself  felt,  in  contemplating  the  exquisite  forms  that  the  camera 
has  revealed,  and  in  watching  the  progress  of  a  multitude  of  events,  compressed  indeed 
within  the  limits  of  a  few  hundredths  of  a  second,  but  none  the  less  orderly  and 
inevitable,  and  of  which  the  sequence  is  in  part  easy  to  anticipate  and  understand, 
while  in  part  it  taxes  the  highest  mathematical  powers  to  elucidate. 

Thus  begins  the  book  A  Study  of  Splashes  (1908)  in  which  A.  M. 
Worthington  (1852-1916)  presents  “in  a  form  acceptable  to  the  general 
reader  the  outcome  of  an  inquiry  conducted  by  the  aid  of  instantaneous 
photography,  which  was  begun  about  fourteen  years  ago.”  Worthington’s 
fascination  with  these  phenomena  actually  went  all  the  way  back  to  1875 
(he  was  then  23  years  old),  when  H.  F.  Newall,  a  student  at  the  famous 
Rugby  boys’  public  school,  gave  a  report  at  the  Rugby  Natural  History 
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Abstract 

A  descriptive  treatment  of  light  emission  from  pulsating  bubbles  (sonoluminescence), 
sound  generation  by  rain  falling  on  water  and  breaking  waves,  and  the  propagation  of 
shock  waves  in  bubbly  liquids  is  presented.  The  first  section  contains  a  bnef  survey  of 
some  occurrences  of  bubbles  in  literature  and  in  the  figurative  arts.  Considerations  on 
the  etymology  of  the  word  are  presented  in  the  appendix. 


1.  INTRODUCTION 

Uon4>6\v^  o  auSpoxcK  -  Man  is  a  bubble  -  is  an  ancient  Greek  proverb  that  enjoyed  a 
singular  favor  in  the  western  culture  as  an  expression  of  the  caducity  and  impermanence 
of  human  life.  It  is  used  by  a  number  of  Latin  and  Greek  authors  such  as  Varro  who, 
apologizing  in  the  Preface  to  his  De  Re  Rustica  (36  b.c.)  that  the  work  is  not  as  polished 
as  he  would  like,  says  that  nevertheless,  being  80  years  old,  he  should  go  ahead  with  its 
publication  since,  “ut  dicitur,  si  est  homo  bulla,  eo  magis  senex”  (“if,  as  they  say,  inan 
is  a  bubble,  all  the  more  so  is  an  old  man”).  Petronius  (1st  century  a.d.),  in  a  mocking 
passage  of  the  Satyricon  in  which  he  compares  man  to  inflated  walking  bap  and  flies,  also 
says  “nos  non  pluris  sumus  quam  bullis”  (“we  are  no  more  than  bubbles,  42,  4).  Lucian 
(117-180)  elaborates:  “I’ve  thought  of  a  simile  to  describe  human  life  as  a  whole  . . .  You 
know  the  bubbles  that  rise  to  the  surface  below  a  waterfall  -  those  little  pockets  of  air 
that  combine  to  produce  foam?  . . .  Well,  that’s  what  human  beings  are  like.  They’re 
more  or  less  inflated  pockets  of  air  . . .  but  sooner  or  later  they’re  all  bound  to  go  pop” 
(^Charon  19). 

'The  advent  of  Christianity,  with  its  message  of  hope  and  salvation,  rendered  the  idea 
less  relevant  and  probably  the  only  pre-Renaissance  textual  reference  to  it  is  in  the  Lexicon 
by  the  lOth-llth  century  Byzantine  scholar  Suida  (“as  a  bubble  immediately  disappeps 
when  it  is  broken,  so  does  the  memory  of  the  splendid  and  powerful  upon  their  death”). 

However,  in  the  northern  post-reformation  cultural  climate  of  the  16th  and  17th  cen¬ 
tury,  the  metaphor  regained  its  appeal.  Most  succint  -  and  first  -  is  Erasmus:  “Homo 
bulla”  [man  (is  a)  bubble,  Adagia,  1508,  n.  1990].  Taverner,  in  his  Proverbs  or  ada- 
gies  ...gathered  out  of  the  Chiliades  of  Erasmus  (1539)  echoes  the  conppt,  and  so  does 
sir  Thomas  Elyot  (1545),  Golding  (“When  man  seemeth  to  bee  at  his  best,  he  is  al- 
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This  paper  presents  the  results  obtained  from  a  detailed  study  of  the  sound  field  within  and 
around  a  cylindrical  column  of  bubbles  generated  at  the  center  of  an  experimental  water  tank. 

' .  The  bubbles  were  produced  by  forcing  air  through  a  circular  array  of  hypodermic  needles.  As 
they  separated  from  the  needles  the  “birthing  wails”  produced  were  found  to  excite  the  column 
into  normal  modes  of  oscillation  whose  spatial  pressure-amplitude  distribution  could  be  tracked 
in  the  vertical  and  horizontal  directions.  The  frequencies  of  vibration  were  predicted  from 
theoretical  calculations  based  on  a  collective  oscillation  model  and  showed  close  agreement  with 
the  experimentally  measured  values.  On  the  basis  of  a  model  of  the  column  excitation,  absolute 
sound  levels  were  analytically  calculated  with  results  again  in  agreement  with  the  measured 
values.  These  findings  provide  considerable  new  evidence  to  support  the  notion  that  bubble 
plumes  can  be  a  major  source  of  underwater  sound  around  frequencies  of  a  few  hundred  hertz. 
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INTRODUCTION 

Wave  breaking  is  the  most  significant  air-entrainment 
process  occurring  at  the  surface  of  the  ocean  and  gives  rise 
to  a  large  number  of  bubble  clouds  (Thorpe,  1982,  1986; 
Thorpe  and  Hall,  1983;  Monahan  and  Mac  Niocaill,  1986; 
Monahan  and  Lu,  1990).  It  is  known  that  the  number  and 
intensity  of  breaking  waves  are  strongly  dependent  on  the 
wind  speed  above  the  ocean  surface  (see,  e.g.,  Toba  and 
Koga,  1986;  Phillips,  1988;  Wu,  1988),  and  it  has  also  been 
shown  that  a  correlation  exists  between  wind  speed  and  the 
intensity  of  low-frequency  (below  1  kHz)  ambient  sound 
in  the  ocean  (Wenz,  1962;  Piggott,  1964;  Perrone,  1969; 
Kerman,  1984;  Kuperman  and  Ferla,  1985;  Wille  and 
Geyer,  1985;  Carey  and  Wagstaff,  1986;  Kennedy  and 
Goodnow,  1990;  Kewley  etal,  1990;  Kennedy,  1992). 
These  circumstances  have  led  to  the  consideration  of  pro¬ 
cesses  by  which  breaking  waves  may  produce  such  low- 
frequency  noise  (Wilson,  1980;  Kerman,  1984).  Carey  and 
co-workers  (Carey  and  Bradley,  1985;  Carey  and  Brown¬ 
ing,  1988)  and,  independently,  Prosperetti  (1985,  1988a, 
1988b)  have  suggested  that  collective  oscillations  of  the 
bubble  clouds  produced  by  breaking  waves  could  be  re¬ 
sponsible  for  the  low-frequency  emissions.  The  argument 
was  essentially  that,  since  the  bubbles  in  the  cloud  consti¬ 
tute  a  collection  of  coupled  oscillators,  one  would  expect 
the  existence  of  normal  modes  of  oscillation  of  the  cloud 
itself  at  frequencies  far  lower  than  the  frequency  of  oscil¬ 
lation  of  the  individual  bubbles.  This  idea  has  been  ex¬ 
plored  in,  and  supported  by,  a  number  of  subsequent  pub- 
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lications  (Lu  etal,  1990;  Yoon  etal,  1991;  Carey  etal, 
1993;  Prosperetti  etal,  1993;  Koller  and  Shankar,  1993; 
Oguz,  1994). 

In  view  of  the  difficulty  in  gathering  oceanic  field  data 
on  the  role  of  bubble  clouds  in  low-frequency  sound  gen¬ 
eration,  the  conclusions  mentioned  above  rest  mainly  on 
theoretical  analyses  only  partially  verified  in  laboratory  ex¬ 
periments.  It  is  therefore  important  to  validate  further  the 
theoretical  models  used  so  as  to  gain  confidence  in  their 
predictions.  Initial  experiments  (Yoon  etal.,  1991)  have 
shown  that  bubble  clouds  are  capable  of  collective  oscilla¬ 
tions  at  frequencies  far  below  those  of  the  individual  con¬ 
stituent  bubbles  in  excellent  agreement  with  theory.  In  this 
paper  that  work  is  extended  in  two  significant  ways.  In  the 
first  place,  the  use  of  a  much  more  extensive  data  set  ren¬ 
ders  the  measurement  of  the  higher  mode  frequencies  pos¬ 
sible  with  a  very  good  match  with  theory.  Second,  the 
model  is  extended  to  the  prediction  of  the  absolute  acoustic 
levels,  again  in  good  agreement  with  data.  This  is  a  very 
nontrivial  point  as  it  presupposes  a  quantitative  under¬ 
standing  of  the  mechanism  by  which  the  bubble  cloud  is 
excited..  Our  conclusion  is  that  the  energy  imparted  to  the 
individual  bubbles  upon  their  formation  coupled  with  the 
spectral  width  of  single-bubble  free  oscillations  accounts 
for  the  level  of  acoustic  radiation  observed  in  the  experi¬ 
ment.  In  a  separate  study  (Oguz,  1994),  it  is  shown  that, 
on  the  same  basis,  good  predictions  of  low-frequency  oce¬ 
anic  ambient  noise  can  be  found.  These  successful  esti¬ 
mates  of  levels  are  a  stringent  test  of  the  theory  as  it  is  well 
known  that,  in  general,  it  is  much  easier  to  match  frequen¬ 
cies  than  levels. 

I.  EXPERIMENTAL  PROCEDURE 

The  experimental  arrangement  used  in  this  work  is 
very  similar  to  the  one  described  in  Yoon  et  al  ( 1991.)  and 
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Bubble  dynamics: 

Some  things  we  did  not  know  10  years  ago 
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Abstract.  This  paper  contains  a  brief  review  of  the  author's  work  on  several  problem  related 
to  bubble  dynamics.  A  “state-of-the-art”  model  for  spherical  dynamics  is  formulated  first.  On 
the  basis  of  this  model  several  features  of  the  thermo-fluid  mechanical  behavior  of  gas  bubbles 
are  then  discussed  and  applications  to  sonochemistry  and  mass  transfer  described.  The  paper 
ends  with  a  brief  discussion  of  pressure  wave  propagation  in  bubbly  liquid  and  the  role  played  by 
bubbles  in  the  generation  of  oceanic  ambient  noise. 

Key  words:  Bubbles,  Bubbly  liquids,  Acoustic  cavitation,  Oceanic  noise 


1  Introduction 

The  period  since  the  1981  lU TAM  Symposium  on  Mechanics  and  Physics  of  Bub¬ 
bles  in  Liquids  [1]  held  at  the  California  Institute  of  Technology  in  Pasadena  has 
seen  a  tremendous  growth  in  bubble  research,  both  fundamental  and  applied.  I 
shall  make  no  attempt  to  review  all  this  new  material,  both  in  view  of  its  amount 
and  of  the  fact  that  most  of  its  originators  are  contributors  to  this  volume.  Rather, 
I  shall  focus  on  some  of  the  results  that  my  colleagues  and  I  have  obtained  in  the 
last  decade. 

I  am  saddened  by  the  fact  that  Prof.  Milton  S.  Plesset,  the  gracious  host  of 
the  Pasadena  meeting,  is  no  longer  with  us.  With  him  we  have  lost  one  of  the 
propulsive  forces  of  our  community,  a  scientist  whose  contributions  have  shaped 
bubble  research  in  our  time.  Personally,  I  am  even  more  struck  as  I  sorely  miss 
his  enlightening  mentorship  and  guiding  voice. 

2  Spherical  Bubble  Dynamics 

Shortly  after  the  Pasadena  meeting  Prof.  Larry  Crum  got  me  involved  in  his  work 
on  bubble  levitation  [2].  In  the  attempt  to  reconcile  the  data  with  theory  we  met 
some  significant  differences  that,  I  am  sorry  to  say,  have  not  been  resolved  to  this 
date.  Around  the  same  time  I  had  also  started  to  work  on  a  book  on  bubble 
dynamics.  Both  factors  prompted  me  to  look  more  deeply  into  the  theoretical 
framework,  and  especially  at  the  incorporation  of  liquid  compressibility  into  the 
well-known  Rayleigh-Plesset  equation,  and  at  the  description  of  the  gas  thermo¬ 
fluid  mechanics.  As  it  turned  out,  answers  to  my  qiiestions  already  existed,  but 
were  not  widely  appreciated  at  the  time. 

As  for  liquid  compressibility,  the  problem  was  not  the  absence  of  equations 
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This  paper  describes  a  theoretical  model  that  predicts  the  wind-dependent  ambient  noise  level  in  , 
the  ocean.  Wave  breaking  and  subsequent  formation  of  whitecaps  are  assumed  to  be  the  sole 
source  of  sound  at  the  sea  surface  and  their  contributions  are  computed  by  the  use  of  a  simple 
model  for  the  bubble  cloud  generated  by  this  process.  Inverted  hemispherical  shapes  for  which 
an  analytical  solution  to  the  wave  equation  is  given  are  employed  to  describe  the  cloud  geometry. 
The  input  physical  parameters  to  the  model  are  the  bubble-size  distribution,  the  dipole  strength 
of  the  entrained  bubbles,  the  cloud  size  distribution  and  growth  rate,  and  the  void  fraction  of  the, 
bubble  cloud.  By  using  an  empirical  relation  between  the  whitecap  coverage  ratio  and  the  wind 
speed,  the  underwater  ambient  noise  and  surface  source  levels  are  computed  as  a  function  of 
frequency  and  wind  speed.  Calculated  noise  levels  are  in  good  agreement  with  the  field 
measurements.  ;  . . 
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INTRODUCTION 

,  The  correlation  between  underwater  ambient  noise 
level  and  the  wind  speed  is  well  documented  in  the  litera¬ 
ture.  (Knudsen  et  al.,  1948;  Wenz,  1962).  In  the  low  (50- 
1500  Hz) -frequency  range,  wind  turbulence,  ocean  spray 
impact,  bubble  plume  oscillations,  and  bubble  cavitation 
are  cited  as  possible  wind-dependent  sources  in  the  litera¬ 
ture  (Kewley  et  al,  1990).  The  cavitation  noise  has  been 
dismissed  as  a  potential  contributor  by  Prosperetti  and  Lu 
(1988)  because  it  requires  unrealistically  large  pressure 
fluctuations  in  the  ocean.  Single  bubble  oscillations  that 
contribute  a  great  deal  to  the  noise  levels  at  higher  frequen¬ 
cies  become  less  of  a  factor  in  this  range  because  it  is 
unlikely  to  find  a  great  many  number  of  bubbles  with  a 
natural  frequency  in  the  order  of  few  hundred  Hz  in  the 
ocean.  However,  the  collective  oscillations  of  bubbles  sug¬ 
gested  independently  by  Carey  and  Bradley  (1985)  and 
Prosperetti  ( 1985)  appear  to  be  a  potential  source  of  noise. 
Experimental  studies  that  have  become  quite  sophisticated 
in  the  past  few  years  all  point  to  the  fact  that  wind  noise  is 
bubble  related. 

: .  Earlier  measurements  of  the  ambient  noise  levels  are 
found  to  be  well  correlated  with  the  wind  speed  rather  than 
the  wave  height  (Penhallow  and  Dietz,  1964)  and  show  a 
logarithmic  dependence  with  a  correlation  coefiicient  that 
is  a;  function  of  frequency  (Piggott,  1964;  Crouch  and 
Burt,  1972).  Wilson  ( 1980)  pointed  out  that  the  noise  level 
follows  the  whitecap  index  with  increasing  wind  speed  at  a 
fixed  frequency  and  attributed  this  behavior  to  spray  im¬ 
pact  noise  (Franz,  1959).  Although  sinlultaneous  ambient 
noise  and  wind  speed  measurements  provide  valuable  data, 
they,  could  be,  contaminated  by  noise  not  related  to  wind, 
especially  at  low  frequencies  where  propagation  loss  is 
small  and  bottom  reflection  high. 

■  ;,In  the  past  decade,  an  innovative  approach  has  been 
taken  to  eliminate  contributions  from  distant 'sources  ,  by 
using  a  vertical  array  of  hydrophones  instead  of  just  one 
omnidirectional  hydrophone.  With  this  technique,  the  sur¬ 


face  source  level  just  above  the  array  can  be  measured 
(Wales  and  Diachok,  1981;  Burgess  and  Kejvley,  1983; 
Carey  et  al,  1990;  Kennedy  and  Goodnov,'  1990; 
Kennedy,  1992;  Chapman  and  Cornish,  1993).  Another 
novel  idea  is  to- visually  record  the  sea  surface  and  the 
acoustic  signal  simultaneously  on  a  video  tape  so  that  a 
direct  relation  between  the  noise  and  the  breaking  waves 
can  be  established  (Hollett,  1994). 

The  collective  results  of  the  more  recent  measure¬ 
ments,  as  well  as  the  old  ones,  make  it  more  clear  that  wave 
breaking,  leading  to  whitecap  formation  and  bubble  en¬ 
trainment,  can  be  a  major  contributor  as  a  sea  surface 
source  in  the  frequency  range  from  50  to  1500  Hz.  The 
modeling  of  the  sound  generated  by  breaking  waves  has 
been  recently  carried  out  by  Loewen  and  Melville  (1991) 
who  considered  only  emissions  from  single  bubbles.  How¬ 
ever,  it  is  well  known  that  the  acoustic  radiation  from 
newly  formed  bubbles  can  be  aSected  by  the  passive  scat¬ 
tering  of  previously  entrained  ■  bubbles  (Updegraff  and 
Anderson,  1991).  Therefore  their  theory  is  limited  to  low 
wind  speeds  when  whitecaps  do  not  exist  in  significant 
numbers.  - 

•  Recently,  a  great  deal  of  attention  has  been  paid  to' the 
understanding  and  measurement  of  whitecaps  and  their 
importance  for  noise  generation.  Detailed  and  careful  ex¬ 
periments  have  been  carried  out  in  the  ocean  to  precisely 
assess  the  role  of  wave  breaking.  Farmer  and  Vagle  ( 1989) 
reported  that  under  identical  wind  conditions  individual 
wave  breaking  events  produce  similar  spectral  noise  levels 
suggesting  a  possible  universality  of  this  phenomenon.  A 
very  detailed  description  has.  been  given  by- Monahan  and 
Lu  (1990)  who  identified  three  types, of  whitecaps  labeled 
as  a  or.  A,  $  ot  B,  and  y.  The  type  A  whitecap  represents 
the  initial  stage  or  the  birth  of  a  whitecap  resulting  from  a 
breaking  wave.  .  It  turns  into  type  B  at  the  end  of  its  life 
which  lasts  about  1-2  s.  Type  A  is  associated  with  the  air 
entrainment  process  and  produces  a  bubble  cloud  that  can 
penetrate  up  to  one  meter  below  the  ocean  surface.  Most  of 
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Abstract.  A  summary  of  the  current  status  of  the  modeling  of  the  propagation  of  pressure  waves 
in  bubbly  liquid  is  presented.  An  apparent  failure  of  the  theory  near  and  above  the  resonance 
frequency  of  the  bubbles  is  found  by  comparison  with  experimental  data.' :  In  the  second  part 
of  the  paper  applications  to  oceanic  noise,  laboratory  bubble  clouds,  and, resonance-parametric 
generation  of  low  frequency  underwater  sound  are  described. 
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1  Foldy’s  theory 


The  lowest-order  theory  of  wave  propagation  in  .bubbly  liquid  was  essentially 
worked  out  by  Foldy  in  a  well-known  1945  paper  [3],  and  later  applied  by  him 
and  Carstensen  to  the  analysis  of  data  [4].  The  following,  although  by  .no  means 
rigorous,  is  a  very  simple  derivation  of  Foldy^s  result.  =  -  -  .  •  ^  '  ' 

In  the  context  of  potential  flow  theory,  to  leading  order,  bubbles  behave  as 
monopoles  with  strength  given  by  r,  the  time  derivative  of  their  volume.  Consider 
then  N  bubbles  immersed  in  an  ‘incident”  flow  given  by  a  potential  <pcc*  The  total 
flow  due  to  the  incident  flow  plus  the  effect  of  the  bubbles  has  then  the  potential 


<f>  = 


N 


<^co  +  X! 

;=i 


47r|x  -  Xjl  \ 


where  Xj  is  the  position  of  the  j-th  bubble.  If  N  is  large,  the  preceding  relation 
may  be  approximated  as 


0  o. 


(2) 


where  n  is  the  bubble  number  density.  Upon  taking  the  Laplacian  of  this  expres¬ 
sion,  since  4>^  is  regular  in  the  region  of  interest,  we  have 

vV  =  nv,  (3) 


which  is  essentially  Foldy’s  result  for  an  incompressible  liquid.  This  can  be  written 
in  a  somewhat  different  form  by  noting  that,  for  the  linear  problem, 

u  =  V(?i>  P  = 

55 
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